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04 ' Using a time-averaging technique we obtain exactly the probability distribution for position and 

^ I velocity of a Brownian particle under the influence of two heat baths at different temperatures. 

^ , These baths are expressed by a white noise term, representing the fast dynamics, and a colored 

• noise term, representing the slow dynamics. Our exact solution scheme accounts for inertial effects, 

' that are not present in approaches that assume the Brownian particle in the over -damped limit. 

, We are also able to obtain the contribution associated with the fast noise that are usually neglected 
by other approaches. 



PACS numbers: 02.50.-r,05.40-a,05.40. Jc 



I. INTRODUCTION 



■ The role of time-scales is crucial for understanding the validity of equilibrium techniques to any physical process 



since fast degrees of freedom (time-scale rj) generally reach their stationary distribution well within the experiment's 
duration but the same cannot be assumed about slow ones (time-scale r^). For instance, glasses are a very important 
class of systems presenting slow dynamics [Ij . 
, The probability distribution for the slow degrees of freedom depends on the system's preparation and also on the 
^ ■ elapsed time interval, i.e., for how long the system has aged For the time interval r/ <C t <C Ts, the system 
O . behaves as driven by an effective external slow field associated with the slow variables. However, there is no physical 
, reason for the slow and fast fluctuations distribution to match, so, in general, Tsiow 7^ Tjast- In this case, the 

' reduced form for the fast variable's distribution is an instantaneous equilibrium one 0- For longer times, t^Ts, the 
T-H I distributions for all variables (fast and slow) reach a stationary state. For closed systems governed by a microscopic 
^ ' Hamiltonian, the stationary state will be a true equilibrium one. For open systems in contact with a single thermal 
On , reservoir at temperature T, the final state will also be the equilibrium distribution. However, by subjecting a system 
■ to simultaneous contacts with two (or more) reservoirs at distinct temperatures, say Ti and T2, it will reach a 
, stationary state distinct from the Boltzmann-Gibbs (BG) equilibrium [s!]. In consequence, the stationary state will 
• ■ depend on the properties of system-reservoir interactions and the reservoirs temperatures. Simple models presenting 
' these characteristics jSuch as Brownian Particles (BP), have been used to explore the physics of slow dynamics and 

: glasses i a i a i i E [o, E [H- 

' The complete description of BP is given by its microscopic interactions with a heat bath of lighter particles as well 
• • as the interactions with the external environment. Trying to evaluate the macroscopic behavior starting up from a 
. ^ microscopic model presents an impossible task. However, for many physical systems we can often reduce the number 
^> , of important variables down to a manageable set, thanks to well separated time-scales [13, [IS EB] allowing us to 
5^ ■ eliminate the fast variables of the problem. The effect of the eliminated variables is taken into account by means of a 
[ random forcing term on the equations of motion for the remaining variables 0, [l7] and by suitable friction coefficients. 
The random term (noise function) represents the effect of the collisions of the BP with the heat baths's particles. 
In the case these thermal forces vary on a very short time-scale, their effect can be represented by uncorrelated 
collisions (white noise) and the BP motion can be analyzed in a simple way because its stochastic dynamics can 
be treated as a Markovian process. Finally, its long term behavior is well described by a Boltzmann-Gibbs (BG) 
equilibrium [Tsj . However, for many realistic Brownian processes in dense fluids the short time-scale approximation 
for the random forces may not be accurate, since the time- intervals of the microscopic collisions might overlap. Instead, 
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we need to consider that those forces act upon the BP during a well defined time-scale, or colored noise. Also, the 
presence of excited large-scale hydrodynamic modes will affect the BP by long time-scale feed-back loops [l^, [13] that 
can be accounted for by a dissipative memory function term in the Langevin-like equation describing the Brownian 
motion (Tsj . 

Our present goal is to obtain the long-time stationary state for a BP system, subjected to contacts with two different 
heat baths at distinct temperatures, that is simple en oug h to be treated exactly but sophisticated enough to present 
the non-trivial stationary properties discussed above j2ll. [2^. 



In our model, it is possible to keep the BP from reaching the BG equilibrium by simultaneously subjecting it to 
the influence of two thermal baths at different temperatures Ti and T2 , one acts as heat source while the other acts 
as a heat sink. As for the BP simultaneously subjected to Ti and T2, the stationary state will depend on the details 
of the interaction between the baths and the BP and the effective temperature will be an intermediary one between 
Ti and T2, as will be shown. Even if temperature gradients are not present, the stationary distribution for velocities 
and positions will differ from the BG measure, reducing to it only when the temperatures are the same, Ti = T2 = T. 
This simple model has been proposed to describe the behavior of glasses [1, H, d, 0, H, M, 113, [iil,[i2,[il[2a,[2l 
where one assumes different time-scales for each thermal bath, via their noise properties: the short time-scale white 
noise reproduces thermal relaxation whereas the long time-scale noise reproduces very slow, correlated, structural 
rearrangement of the glass. Additionally, in our approach we can include the effects related to the finite mass of the 
particle in stationary state, as well as other effects associated with the slow noise Q which are often neglected by 
other methods based on the over-damped limit calculations (which corresponds to the limit m —^ 0). For instance, 
the BP adjusts itself to the slow noise function as [3, HI], specially when its mass can be neglected (over-damped 
case), and the position distribution becomes very similar to a true equilibrium form, given that the slow noise acts 
like an external effective potential Q . However, this simplification eliminates the contribution of the fast noise for 
the stationary distribution fS]. In contrast, our time- averaging scheme will integrate that small effect for a very long 
time and exhibit the missing contribution. 

This paper is organized as follows. In Section II we describe the model. In Section III we calculate the time- 
averaging integrals obtaining the displacement and velocity non-vanishing terms. In Sections IV and V, we obtain 
the stationary (or equilibrium) distribution and discuss our main results. 



Let's consider a massive Brownian Particle (BP) moving under the action of a confining potential V{x), which is in 
contact with two heat baths of distinct temperatures and time-scales, expressed by the noise functions ^(t) and ri{t). 
There are some interesting physical systems that present this kind of multiple time-scale phenomena. For instance, in 
glasses the short time-scale (fast) modes are associated with thermal vibrations and the long time-scale (slow) ones 
with structural reorganization of the atomic structure [2^ . 

The effect of the fast modes on the BP are represented by the white noise term 77 (i), associated with the heat bath 
at temperature Ti and the instantaneous friction coefficient Pi. On the other hand, the slow modes contribution are 
represented by a colored noise term £,{t), associated with the heat bath at temperature T2 (in general T2 ^ Ti). These 
modes give rise to long-time feed-back loops that are expressed by the integration of a frictional memory function 
(/)(t — t ) over time. 

This is modeled by two coupled stochastic differential equations: 




II. EXACTLY SOLVABLE MODEL 



A. Langevin-like Equation 



X — V, 



(1) 




(2) 



where the properties of the terms above are discussed in the following. We assume V{x) 



^kx^ for simplicity. 



B. Noise properties 



The stochastic process described by Eq.© is non-Mar kovian due to the noise ^ be time correlated with a well 
defined time-scale r and exponential time-correlation function, or a colored noise (see for example ,25.]*). It represents 
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the slow dynamics and only at time-scales much larger than r will ^ manifest Markovian behavior. The noise 77 is a 
white one and represents the fast dynamics. Both noises, fast and slow, are Gaussian and can be defined in terms of 
their two lowest cumulants. Let's start with the slow noise properties: 

m) = 0, (3) 

D2 ( \t~t\\ T2T2 ( \t-t\\ 

mat )) = ^ exp = ^ exp , (4) 



giving the memory function ^26*1 



,s mat )) r2 / \t-t\ \ 

(j){t-t) = — = — cxp I — 1 , (5) 



where t is the characteristic time-scale of the slow noise, T2 its temperature and T2 a dissipative strength associated 
with ^. The units are chosen so that Boltzmann's constant is equal to one (fcs = 1). 
The first two cumulants for the fast noise are given below: 

m) = 0, (6) 

{r^{t)f^{t))=Tir,5{t^t'), (7) 
where Ti is the temperature and Fi is the dissipative coefficient for the fast noise. 



C. Equilibrium distribution 

Our goal is to find out the exact stationary solution for Eq.Q through the time-averaging of 

p{x, V, t) = 5{x - x{xo,Vo, t))5{v - v{xo, Vo, t)), 

as t —* 00. The functions x{xo, Vo, t)) and v{xo, Vq, t) are the solutions for a given set of BP's initial conditions {xo, Vq) 
and for a given realization of the stochastic processes (,{t) and ri{t) We need to take the average over these last 
two stochastic processes in order to obtain the stationary solution. 

Thus, the stationary distribution for the Brownian degrees of freedom reads [2^ : 

1 

P'%x,v) = lim - / dt {p{x,v,t)}^ 

poo 

= lim z / dte-'' {p{x,v,t)), (8) 

z^0+ Jq 

For simplicity, we will assume that the initial conditions are Xo — Vo ^ 0, since the terms depending on the initial 
conditions tend to zero as t ^ 00 [2611. 



D. Laplace transformations 

The centerpiece of our method is to use Laplace transforms in order to make our problem manageable and the 
mathematical reason for doing so is because Eq.® itself is the Laplace transformation of the averaged density. Thus, 
our goal is to describe the stationary state probability distribution, P''''{x,v), in terms of the Laplace transform of 
the noise functions, |^(z) and fj{z) [26j . 

In fact, we can express P^^{x,v) as a sum of many terms involving ^{z) and fi{z). Most of these terms will not 
contribute to the distribution, but we can calculate analytically the remaining ones, finding the exact expression for 
P^''{x,v). In order to do that we need to express the Laplace transforms, of the position and velocity of the BP, as 
functions of ^(z) and 77(2:). In the following we proceed to implement the strategy just described. 
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We start by taking the Laplace transformations of Eqs.JT]) and ([2]) (with Re(z) > 0): 



Ta k 
mz ^ + Fi H — 

T Z + 1 Z 



and 



zx(z) = v{z). 

The Laplace transforms for the noise correlations are given by [33j : 



(9) 



(10) 



+ e) iiiqj + e)) 



2r2T2 



i{q, + qj) + 2(] [1 - T{iq^ + e)] [1 - T(ig, + e)] 

2, + T[i{q, + qj) + 2e] - T^{iq, + e){iqj + e) 



T2T2T 



[1 - T{tq, + e)] [1 + T(zg, + e)] [1 - r(ig, + e)] [1 + T{tqj + e)] 



(11) 



2riri 



i(qi + qj) + 2e' 



(12) 



(??(«g« + e)^(«g« + e)) = 0. 



(13) 



It can be shown that the second member of the RHS of Eq. (fTT|) does not contribute to the stationary distribution 
P^''(x, v) when z ^ in Eq.®. The above Eqs.Hl]), ^ and ^ will be crucial for calculating P'"*(x, v) via Eq.®. 

We observe that when Fa — *■ we recover previous results (26l ) since the contribution of the slow noise disappears. 
Furthermore, we notice that by taking the limit t —^ the colored noise ^ becomes a white one, which is reflected in 
the fact that its contribution in Eqs.([2]) and ([9]) becomes similar to the one by rj. 

Following a strategy similar to the used in Ref. ^2^ , we express the Laplace transformation of position and velocity 
as functions of the Laplace transformation of the noise functions: 



v{z)^e{z) + , 



(14) 



with 



where 9 and u are given by 



and 



e(z) zn{z) = 



i{z)^n{z) [az) + v{z) 

z{tz + 1) 



V mr niT rar I 



TO + Fi T, 



= Fi + F2 + A: r. 



(15) 



(16) 



(17) 



(18) 



For simplicity, we assumed the initial conditions xi^ = 0) = and v(t = 0) = 0, since for T 00 the terms carrying 
the memory of the initial conditions will decay to zero [26[. The denominator of 6(2;) and r2(z) is of the third order 
on z, thus, there are three distinct roots for: 



z-" H z'' + z + = 0, 

TOT TOT TOT 



namely zi, Z2 and Z3: 



Zl 



3 TO T A 



(19) 



(20) 
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22 = ^3 = -^ — ^ — ^ — +'—F — r^' (21) 

brriT A brriT \ 

where 

a= + 3ujmT, (22) 

and 

mr [3 (4 w^mr — — 18 uj 9 rriT k + 27 k^rriT^ + 4 fc6'3)]i/2. (23) 

We define the quantities A, B, C, V, T and Q as: 

A^iB \ z^^ z-2. = C^iV; z^^ z-2. = T ^iQ. (24) 

The terms A, B, C, V, T and Q are important because we perform residue integrations over their poles [l^. It is 
important to compute the expressions for the leading terms in the expressions ((22)) . ((23)l and ((24|) when mr ^ 0, i.e., 
in the over- damped limit. Since > 0, in linear order on mr we have 

, / iLoniT mT\/3 I \ 

A = {-9) h - — _^ v^^^nnj . (25) 

The function o is already linear in m t. In the leading order in m r, we obtain: 



^ TT^' ^ ' • 

29 20 mr 20 niT 20 

E. Calculating the stationary distribution 

In Appendix |A] we show the derivation for the following equation (originally obtained in Ref. [26]) for the stationary 
state distribution Eq. (S) 

P-(x v) - hm V '^Pe^Q^+^P^ ^-'^^ ("^-P)"' (^^ A '^'^^ TT '^^^ 

[^x,v) - lim 2. 2^ e y YV 27r -li 27r 



n ^(*'?/ + ^) n + ^)) (27) 



In order to compute P'** we need to replace Eqs. ([TT|) to (fT5|) . into Eq. ([27|) . The integrations paths for the variables 
{qf,Ph} can be chosen to be identical to the one shown in Fig.([T]): the integrations in Eq. ((27|) will correspond to a 
series of residue integrations over (only some) of the multiple poles of its integrand. These poles correspond to the 
ones associated with the noise functions fEas. (fTT| ) to (|15p ) and the ones given by Ea. p6)l and the function I{z): 

Hz) = — -■ (28) 

E f=i "ilf + EhLi "^Ph + + m)e 

In special, the function I(z) is crucial for our computations but its role only becomes clear near the end of the 
calculations. As a matter of fact, the pole of I(z) is located outside of the integration path for all {qf,ph}, as in Fig.([T]). 
The really important poles are associated with the roots of the denominators of Q{z) and Q{z), zi, zi and Z3 above 
(the z-poles), and the ones from the denominators of the Laplace transforms for the noise variances (the noise-poles), 
Eqs.fTTjl and (fT2|) . that are locate on the inside of the path of Fig. ([1]). After replacing the rather complicated forms for 
the noise Laplace transforms into Eq. ([27| . the computing of P'**(x, ti) becomes a task of doing the residue integrations 
around the relevant poles and collecting the non-vanishing terms. 

The effect on l(z) due to the residue integration of Eq. ([27|) over its relevant poles is rather interesting. The 
integrations over the noise-poles for that are present in l(z) imply that the value the denominator of I(z) 

will take has a non-zero real part, even if only one of such an integration has been done. In that case, when taking the 
limit z ^ 0, we obtain limz^Ql{z) = 0. However, whenever an integration over any one of the {qf,ph} variables at 
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a noise-pole is done, it has the effect of reducing a sum of the form i{q+p) + 2e, i{q + q') + 2e or i{p-\-p') + 2e to zero 
in the denominator of I{z). It means that the other member of the pair may be integrated around any of its poles 
without contributing to the denominator of /(z), making that integral a possible non-zero contribution to P'"'{x,v), 
when the limit z — > is taken. Thus, we notice that only the residue integrations that eliminate the {qf,ph} variables 
from the denominator of /(z), two of them at a time, will reduce I{z) to 

2; 

limz^al{z) = liniz^Q— = 1. 

z 

Only in this case, I{z) — 1, the corresponding term will contribute to the stationary state distribution 26]. 

It is straightforward to show that cross terms that couple x and v (il0-type) do not contribute, while other terms 
coupling X and x (il57-type) or v and v (00-type) will give non-zero contributions. As a consequence, the stationary 
probabilities for positions and velocities are independent: — P'"*(a;)P'*'*(t;). 

In the following, we evaluate some typical integrals and deduce the stationary distribution. 



III. TIME-AVERAGING 



A. Inertial case 



In this section, we evaluate the terms that generate the stationary state distribution P'^'^{x,v) by analyzing 
the distinct contributions grouped into Eq. (P7)l . So, lets start by studying a pair- integration contribution from a 
rifi— correlation integral (for a position-position correlation this is a [q, q'] pair, such as mentioned in the last section): 



dqi dqj 



27r 27r z — i(qi + qj — 2ie) 



n{iqi + e)n{iqj+e) (fi{iq, + e)fj{iqj + e)) + + e)f(igj -|- e)) 



(29) 

where we integrate over the contributing poles (only), and io — ^ j iqf + + {I + m — 2)e represents 

the summation over the p's and g's that are not integrated over in the denominator of I{z). The Wn factor is given 
below (see Appendix [C]) : 



TiTi 2ziAB+{BVg -BCT) + {AVT + ACg) 
2m2 y^S(C2+I?2)(jr2 + g2) + 

{T^T^+T2T2) 2AB{A^ +B^)- z^{BVg -BCT)+z^{AVT + ACg) 



2 m? t2 



Zi AB{A^ + 62) (C2 + I?2) (_^2 + ^2) 



(30) 



The terms depending on Fi Ti (r2 r2) result from the integration over the poles of the {fjf}') and the effect 

of the integration over qi and qj on Eq. (^5]) corresponds to the appearance of a Wn factor and the reduction of the 
denominator of /(z) from [z — i{qi + qj — 2ie) — io] to [z — io]. 

Similarly, a typical pair-integration contribution from a 00— correlation integral (velocity-velocity correlation — > 
pair) reads: 



dpi dp. 



— Q{ipi + e)e{ipj + e) {f]{ip^ + e)fj{ip.j + e)) + (^(ip^ + e)S,{ipj + e)> 



2ti 271 z — i[pi + pj — 2ie) — io 
where the Wq factor is 

(FiTi +r2r2) 2ziAB+{BVg -BCT) + {AVT + ACg) 



z — o 



-We 

> 

(31) 



We - 



2m2r2 AB{C^ +V^){T^ + g^) 

TiTi 2zfAB + {A^ -3AB^){Cg + + {B^ - i A^ B){C T ~ V g) 



2m? 



AB{C^+V2){T2 + g^) 



(32) 



It is straightforward to show that the integrals of the type Ofi— correlation (for a [q,p] pair) over the correct poles 
will not contribute to P*^(x,u). It means that the information about the position distribution (g's) does not get 
mixed up with the information about the velocity distribution (p's). In other words, the distributions are independent 

(P*''(x,ti) ~ P^'*(a;)P'**(w)) as mentioned before. 
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Schematically, we can represent the integral in Eg. ((27)) by: 

2n 2m all p-pairs all g-pairs , . 

linrJn^^.n^'^'T^ n fe^.) n fex.)->(jnrs)^e^- 
j=i 1=1 ij ki ^ ^ 

where x represents the noise functions' Laplace transforms. The number of terms above remains to be evaluated. 

First we notice that we can factor out the Wq contributions from the Wq contributions. By doing so, we only 
need to calculate the number of ways of integrating all [q,q'] pairs (and similarly for the [p,p'] pairs). Furthermore, 
averaging over a power 2t of ^'s and 77's yields a number of terms analogous to the number of ways of distributing 2t 
balls into t boxes that have room for two balls each. The factor above is given by: 

(2t)! 

Let's replace that into the equation for P'"*(a;). The factored out contribution for the position distribution reads 

°° /•+00 1^ f nn\'^i 



h'-oo 2- (20! 2'/!^"' 



P-(a;) = J_^ e-^. (33) 
V2 7r Wn 



Similarly, it is straightforward to show that: 

P^^iv) = J_^ e-^. (34) 

B. Stationary distribution 

The Eqs.(I2Sl) and ^ combine to 

P"'(x,v) = P'''(x)P'''{v) = / ^^pj-- (— + — n . (35) 

The distribution on Eq.(|35p is complete, normalized and exact. It represents the stationary state for extremely long 
times and includes more information than similar models in the massless (over - damped) limit. Since Wn 7^ Wq, a 
renormalization of the mass and of the t emp erature will occur. Recent works have suggested other mechanisms for 
deviations from Boltzmann statistics [l^, l28l . [29j . A massless BP is bound to move as the instantaneous force directs 
it, but an inertial one can accumulate kinetic energy and its velocity deviates from the ratio between the external 
force and the friction coefficient, creating an opportunity for a feed-back mechanism (incipiently present in the non- 
Markovian colored noise ^ time-scale r) to modify the stationary distribution away from the Boltzmann equilibrium 
form [2^, [33, [31I, [12] even at Ti = T2 pg . [soj . Taking the limit ni ^ at Ti — T2 restores the Boltzmann equilibrium 
form, despite r > 0, due to the crucial role of inertia in the mechanism above. 

IV. OVER -DAMPED LIMIT 

The Eg. ([30)1 defining the Wq term, (see Appendix |B] for a complete derivation), gives the dependence of the 
stationary position distribution on all parameters of the system in its most general form. This distribution is a 
generalization of the expression obtained in Ref. [2^ and, similarly to that case, it differs from the Boltzmann form 
due to the non-Markovian character of the colored noise, but tends to it when we take r ^ and Ti = T2. The 
nontrivial forms of Egs. l|5D)) and indicate that a renormalization of the mass (or the rigidity k) takes place due 
to the non-Markovian character of ^ [2^, [Sll, [s^ • 

In the over damped limit, m 0, the term Wn simplifies considerably and reads 

Ti r fc Fi Ti -h Fz T2 
lim Wn^ \^ , 36 

m->o (f 1 + f 2 + kr) k 
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The quantity kWci behaves as the exact effective temperature T^ff for the system in the over -damped hmit. We 
observe that both the slow and the fast noise contribute to T^fj. It can be seen that T^jf — k Wq has an intermediate 
value between Ti and T2: 

min[Ti,T2] < kWn < max[ri,r2]. 

Comparing Ea. ((36| with other approaches, we noticed that in Ref. Q the authors show the steps to obtain the 
stationary distribution for a similar model in the over damped limit. In their case, they make approximations in such 
a way that both the slow field ^ and the dissipative memory term, proportional to r2 (Fi in their notation), behave as 
if they were effective external additive terms for the potential V{x). In consequence, the frictional coefficient Fi (Fo 
in their notation) will not contribute to the equilibrium distribution since that is a function of the effective external 
terms only. The system is supposed to quickly adjust to that potential, giving rise to an instantaneous equilibrium 
distribution. 

On the other hand, in the present method we calculate the averages rather differently. As we have seen, we first 
integrate over the noise-average of the exact distribution, Eq.®, for an infinitely long-time and then obtain the 
time-average of them. By doing this, the neglected effect of Fi becomes apparent and shows in our solution. So, in 
order to recover the results suggested on section 4.1 on Ref. Q (derived in appendix ICj) all we need to do is take 
Fi = in Eg. ([55)1 . In our notation, the results on Ref. '3i] read (where we assume V{x) — so our models coincide) 



W, 



[15] _ TiTk + T2T2 

{T2 + kT)k ' 



n 



which is identical to ours when we take Fi = in Eq. ([36|) . 

When Ti = T2 = T we recover the Boltzmann equilibrium with distribution 



W^o = ^^P-(x) = ;^e-^. (37) 

For the velocity distribution, in the over -damped limit it assumes a simple form. As is the case of Wn, the 
expression for Wq in Eq.([32l) expresses the dependence of the velocity distribution on all parameters of the model. 
However, by taking m — > in Eq. ([32|) gives asymptotically: 

lim We ^ —. (38) 

m— >0 TO 

We notice that only the fast noise Ti contributes to the velocity distribution. The kinetic energy of the BP is 
completely driven by the fast noise and ignores the slow one in the over -damped limit. 
The normalized velocity distribution reads in that approximation order: 



We^ — ^ P'^iv) = /— l^-e"^. (39) 



In this limit we observe that only the fast noise (Ti, Fi) drives the kinetic energy distribution. 



V. CONCLUSIONS 



A better understanding of slow dynamics and associated models may require powerful numerical methods, and 
computers, as one tries to simulate the (extremely) long-time relaxation that often occurs for interesting systems such 
as glasses. 

In this context, simpler models may be very useful for obtaining some good qualitative understanding of the long 
time limit without presenting the numerical difficulties the more realistic ones do, with the possibility of obtaining 
exact results. These are always interesting in that they take into account all the physical effects present in a given 
model, independently of any approximation. In other words, exact treatments contain all the available information 
about a model. Some of it being inaccessible through approximate methods. 

In that spirit, BP models have been proposed as a means to study some of the phenomena associated with the 
competition between fast thermal fluctuations and slow structural relaxation in glasses [1, [2^, A white noise 
(fast) is associated with the thermal fluctuations that happen in short time-scales. A colored noise function (slow) is 
associated with the long time relaxation (structural). 
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For these simple models, a typical approach is to consider the over-damped regime, which is equivalent to the case 
when the BP has zero mass (m — !■ 0+). When responding to external forces acting on it the BP's velocity assumes 
the ratio between the sum of all external forces and a friction coefficient. These models can be used to study the 
emergence of a non-equilibrium stationary state as the probability distribution tends to the stationary form as t ^ oo. 

In the present work, we develop an exact method for time-averaging the distribution for x and v for a BP submitted 
to a fast white noise and a slow non-Markovian one at different temperatures. 

The essence of our method is the use of time-averages of the distributions for x and v which leads us naturally to 
the use of Laplace transformations in order to "open" the problem. In a straightforward way, we represent the time- 
averaged distribution as a sum of integrals that can be easily analyzed, and computed exactly. The distribution thus 
obtained can be compared to the ones in the literature. By being exact, our method allows us to obtain information 
about the effects due to the finite mass of the BP showing, for instance, that the fast noise term plays a role on the 
stationary distribution of positions that is usually neglected in other methods. 
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APPENDIX A 

The Laplace transformation of Eq. ((27|) is shown below [2^ . We will give a step by step derivation for the expression 
of P'"'{x, v) as a function of the Laplace transforms of the position and velocity for the Brownian particle. 
Lets start with the basic definition: 

/■oo 

P'''(x,v) = lim z / e~'US(x ~ x(t))S(v - v(t)))dt. 

Jo 

We write the delta-functions above in the integral representation: 

{Six - xitmv ^ vim = r^e^'^^fl^-^f:^-^ 

and obtain (after using the delta functions to express identically the averages over the noise as functions at distinct 
times): 

P-(x,«) = lim z r dte-^' ^°°^^^Q^ f^°° dP^.P.yt^y 



z^o+ In /—no 27r /—no 27r 

/ YldtifSit-tif) Y[dt„,hSit-t,nb){Y[xitif)Y[v{t^h)). 

•'^ f=l •'^ h=l f=l h=l 



f=l h=l f= 

Next, we express all delta functions above as integrals on the complex plane, displaced from the complex axis by 
a factor of e (that vanishes faster than z) . That factor will guarantee the convergence of the Laplace transforms for 
positions and velocities in the following. 

n+oo r + oc ir, oc , -^w 



-.^-0+ Jo J^oo 27r 2TT ^ n ml 

Jo "'O 



I 



We need to integrate over all {t;/,tmft}j obtaining the averages over the Laplace transforms of the position and 
velocity: 

-^-^ 1=0 m=0 /=1 h=l 

X / d^ze-*{^-^/=l('"'^+')-^"=l(*^"•+")}(]^^(^g/+e)]JD(^p,, + e)). 

•^0 /=! h=l 



Finally, we integrate over t and obtain Eq. ([?7|) : 

l,m=0" 



eLi + Er=i + + ™)'^ I /=i ''^i 
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APPENDIX B 

A typical calculation of the stationary state distribution of displacements terms is shown below for the Wn term 
from Eq.dlSl). 

The first integral in Ea. (P5)) reads: 

TTIT' 7 \ — TTT + e)n{iqj + e)(?7(% + eMiq^ + e)) = 

J ZTT ZTT z — + — 2^ e + o) 

^ dqi dqj z 1 + t (iqi + e) 



27r 27r z - + g^- - 2i e + o) [g, - i(e - zi)] [g^ - i(e - Z2)] - i(e - Z3)] 
1 + T {iqi + e) Ti Ti m-2 



[qj -z(e-zi)] -i(e-Z2)] - i{e - Z3)] {-i)[qj - {-qi + 2ie)] 

To continue the calculation, the integrations over the poles must be done in a way that allows us to obtain, after all 
integrations have been done (as explained in Section III. A), lim2_»o = 1- When integrating over the poles of qj's 
it is possible to see, Fig.dT]), that it will work only for qj ~ —qi + 2ie. Thus: 

^ 2^2^z-z(g. + g,-2ze + o)"(^* + '^^^'''^ + ^^^^^^'^'^ + '^"^^"^^ + = 

FiTi z r+°° dq^ 1 - {iq, + e)^ 

t2 z - io 7_oo 27r [9^ - i(e - zi)] [g^ - i(e - Z2)] - i(e - Z3)] [g^ - i(e + zi)] [g^ - i(e + Z2)] [qi ~ i{e + Z3)] 

The same holds for the integration over the q,;'s poles. So, the non-zero contribution comes only from the poles 
qi — z(e — Za), a — 1,2, 3, in the upper part of Fig.p]): 



dqi dqj z 



■ Q.{iqi + e)fl{iqj + e){fj{iq^ + e)f]{iqj + e)) 



27r 27r z — i((7i + (7j — 2i e + o) ' '^'"'^"^^ * 

FiTi Z r t2z2-1 r2[z2(zi-Z2)*(zi+Z2)*-Z2*(zi-Z2)(zi+Z2) 



z — io [ zi I zi — Z2 p I zi + Z2 p 4 i 3fi(z2) 3(z2) | zi — Z2 p | zi -t- Z2 p 

[z^(zi - Z2)*(zi + Z2)* - Z2(zi - Z2)(zi + Z2)] 



4Z3?(Z2)3(Z2)|Z2P|Z1-Z2P|Z1+Z2P 

Finally, using the definitions of Eq. ([24l) . we get for the first part of Eq. ((29l 
dqi dqj z 



Q.{iqi + e)fl{iqj + e){fi{iqi + e)T){iqj + e)) = 



^ 27: 27: z — i{qi + qj - 2i e + o) 
z (TiTi 2ziAB + {BVg - BCT) + (AVT + ACg) 



Fi Ti 2AB{A^+B^)~zi{BVg-BCJ') + zi{AVJ' + ACg) \ 

For the second part of Eq. , the second term of the RHS of Eq. (fTT|) docs not contribute since it is straightforward 
to show that it leads to a null contribution. Thus, we shall compute the contribution from the first part of Eq. pQ]) 
keeping only the RHS of Eq. (fTT|) . It yields: 

, 271 271 z — i[qi + qj ~ 2i e + o) 

^ dqi dqj z 1 + r [iqi + e) 



27T 27r z - i{qi + qj - 2i e + o) [qi - i{e - zi)] [qi - i(e - Z2)] [qi - i{t - Z3)] 
1 + T {iqi + e) F2 r2 m-2 



- «(e- zi)] [qj - i(e- Z2)] [qj - i{e- Z3)] {-i)[qj - + 2ze)][l - riiqi +e)][l - T(igj- + e)] 
The only contributing integrations will be the ones over the poles qj = —qi + 2ie, (see Fig.([T])). So: 



(Bl) 
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TaTa z dq^ 1 1 



, 27r [qi - i{e - zi)] [qt - i{e - Z2)] [qi - i{e - Z3)] [g^ - i(e + zi)] [qi + i{e - Z2)] [qt + i{e - Z3)] 
Again, the contribution came from the poles qi — i{e — Zq), a = 1, 2, 3. So: 

ZTT 27r z — i[qi + qj — Zi e + o) 

T2T2 Z ( -1 [Z2(Z1 - Z2)*(2:i + Z2)* - Z2(Z1 - Z2)(Z1 + Z2)] 



2m2r2 z- io zi | zi - Z2 p | zi + Z2 p ' 4i 3?(z2) 3(z2) | Z2 p| zi - Z2 p | zi + Z2 p 
Using the definitions on Ea. ([M)) we obtain: 

00 27r 27r z — + — 2i e + o) 

r2r2 2AB{A'^ + B^)-zi {BVg -BCT) + zi {AVT + ACg) 



2m2T2 ziAB{A^ +B-^){C-^ +V^){T'^ + g'^) 



(B2) 



Combining Egs. ljBip and (jB2p we obtain the resuh for Wji in Eg. ([50)1 . An analogous calculation can be done for 
the velocity integrals and obtain We, e.g., Eq. 
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APPENDIX C 

The stationary distribution of x is obtained in Cugliandolo and Kurchan Q. We review their main steps below. 
Let's define: 

P{x) = J dhP{x/h)P{h), (CI) 

where 

P{h) = — (C2) 

Jdhe 

and 



Pixlh) = . (C3) 



The denominator above defines Z{h) and F{h) = —f3 ^ In Z{h). For V{x) = fc^, we have: 



Z[h) = , / e^<'=+^' ; F{h) = ^ p-'C (C4) 



And the denominator of Ea. (jC2[) becomes: 



rf/^e V .) =e ^ W (C5) 



Thus, Eq.fCTj) becomes: 



/ 3* 3 k(kT + r9) f^' ^ fc(fcT+r2)x^ 



When the temperatures are equal, then (3* = f3 = 1/T and: 
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